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Abstract 

Esophageal cancer remains the sixth leading cause of cancer associated death and eighth most common cancer worldwide. 
Genetic factors, such as single nucleotide polymorphisms (SNPs), may contribute to the carcinogenesis of esophageal 
cancer. Here, we conducted a hospital based case-control study to evaluate the genetic susceptibility of functional SNPs on 
the development of esophageal cancer. A total of 629 esophageal squamous cell carcinoma (ESCC) cases and 686 controls 
were enrolled for this study. The OPG rs3102735 T>C, rs2073618 G>C, RANK rs1805034 T>C, RANKL rs9533156 T>C and 
rs2277438 A>G were determined by ligation detection reaction method. Our findings suggested that RANK rsl 805034 T>C 
is associated with the susceptibility of ESCC, which is more evident in male and elder (>63) patients. Our study provides the 
first evidence that functional polymorphisms RANK rsl 805034 T>C may be an indicator for individual susceptibility to ESCC. 
However, further larger studies among different ethnic populations are warranted to verify our conclusion. 
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Introduction 

Despite recent considerable medical advances, esophageal 
cancer remains a refractory disease with high morbidity and 
mortality. Essentially, esophageal cancer is the 6 th leading cause of 
cancer-related mortality and the 8* most common cancer 
worldwide [1]. There are more than 450,000 patients diagnosed 
as esophageal cancer worldwide and the incidence is still rising 
rapidly. Meanwhile, its startling overall 5-year survival rate ranges 
from 15 — 25% [2]. In China, more strikingly, esophageal cancer 
ranks the 5 th most common diagnosed cancer and 4 th leading 
cause of cancer related mortality [3] . Esophageal squamous cell 
carcinoma (ESCC) is the predominant histological type of 
esophageal cancer [1]. Although multidisciplinary therapeutic 
strategy has been recommended, the prognosis is still poor. 
Tobacco use [4,5], alcohol consumption [4,6], low socioeconomic 
status, poor oral hygiene and nutritional deficiencies [2,7-9] have 
been identified as risk factors for esophageal cancer. Yet, only a 
subset of individuals exposed to these risk factors eventually 
develop esophageal cancer, indicating a pivotal role of genetic 
factors, such as single nucleotide polymorphisms (SNPs), in the 
esophageal carcinogenesis. 

Recently, the osteoprotegerin (OPG), its binding protein-the 
receptor activator of NF-kB (RANK) and RANK ligand (RANKL) 
have been implicated with the pathogenesis of breast cancer [10]. 



OPG was initially identified from a fetal rat intestine cDNA library 
[11], which is unique for it only exists as a secreted molecule in 
contrast to the other membrane-bound cell surface members of 
tumor necrosis factor receptor (TNF-R) family. RANKL is the 
OPG binding protein (also named OPG ligand, OPGL) [12,13], 
while RANK constitutes the cell surface receptor which responses 
to OPGL. In numerous rodent models of tumor, RANKL signal is 
increased through diverse mechanisms [14]. OPG neutralizes 
RANKL, which leads to a reduced RANKL-RANK interaction 
[12]. RANKL expression was verified in various tumor types and 
inflammatory cells associated with tumor [15—17]. Elevation in 
stromal RANKL has been detected at local sites of bone metastasis 
or multiple myeloma [18,19], causing enhanced osteoclast activity 
and bone destruction. In experimental models, RANKL inhibitors 
reduced tumor-induced osteolysis in various types of cancer [14], 
reduced bone destruction, skeletal tumor progression, as well as 
tumor burden [17,20,21]. In addition, RANKL-RANK pathway 
may contribute to the primary tumorigenesis and metastasis 
independently of its effects on tumor-related osteolysis. Regulated 
by factors including prolactin and progesterone, RANKL could 
drive the primary mitogenic response of mammary epithelium and 
the expansion of mammary stem cells via RANK activation [22- 
24], which may therefore induce mammary cancer by offering a 
more transformation-susceptible target pool. RANKL may regu- 
late spontaneous mammary tumor formation and metastasis 
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driven by the potent oncogene Neu (ERBB2). RANKL blockade 
effectively attenuated the formation of mammary tumors and 
pulmonary metastasis in the MMTV-Neu transgenic mouse model 
[25,26]. Interestingly, OPG may serve as a positive regulator of 
microvessel formation and may promote neovascularization [27] 
that is important for tumor progression. OPG overexpression by 
breast cancer cells enhances orthotopic and osseous tumor growth 
[28]. In light of all these findings, RANKL/RANK/OPG 
signaling pathway has emerged as a promising therapeutic target 
of cancer. Denosumab, a monoclonal antibody against RANKL, 
has been approved for the treatment of postmenopausal osteopo- 
rosis and bone metastasis in breast cancer [29] . 

Genetic variations in genes encoding RANK, RANKL and 
OPG were found to affect the rheumatoid arthritis [30], Paget's 
disease of bone [31], hip osteoporotic fracture [32]. More 
importandy, in the Caucasian population, a significant association 
of the SNP rs3 102735 (OPG) with the susceptibility to develop 
breast cancer has been reported [10]. The functional significance 
of RANK/RANKL/ OPG signaling pathway and pilot study in 
breast cancer have led us to investigate the association between the 
esophageal cancer and SNPs in the genes oiRANK, RANKL and 
OPG. In a hospital-based case-control study, we performed 
genotyping analyses of the live miRNA SNPs in 629 ESCC cases 
and 686 controls in a Chinese population. 

Methods 

Ethics Statement 

This hospital-based case-control study was approved by the 
Review Board of Jiangsu University (Zhenjiang, China). All 
subjects provided written informed consent to be included in the 
study. We have complied with the World Medical Association 
Declaration of Helsinki regarding ethical conduct of research 
involving human subjects and/ or animals. 

Study populations 

A total of 1 ,3 1 5 participants consisting of 629 esophageal cancer 
patients and 686 non-cancer controls frequency-matched to the 
cases with regard to age (±5 years) and sex were enrolled in this 
study (Table 1). All patients and controls were consecutively 
recruited from the Affiliated People's Hospital of Jiangsu 
University and Affiliated Hospital of Jiangsu University (Zhen- 
jiang, China) between October 2008 and December 2010. All 
cases of esophageal cancer were diagnosed as ESCC histologically. 
Patients who had cancer history/metastasized cancer or had 
received chemotherapy/ radiotherapy were excluded for the 
current study. 

Each subject was personally questioned by experienced 
interviewers using a questionnaire to obtain information on 
demographic data (e.g., age, sex) and related risk factors (including 
tobacco use and alcohol consumption). After the interview, 2-mL 
samples of venous blood were collected from each subject. 
"Smokers" subgroup included individuals who smoked one 
cigarette per day for >1 year. Subjects who consumed ^3 
alcoholic drinks a week for >6 months were subdivided into 
"alcohol drinkers" category. 

Genomic DNA extraction, SNP selection and Genotyping 

Genomic DNA was isolated from peripheral blood using 
QIAamp DNA Blood Mini Kit (Qiagen, Berlin, Germany) as 
reported previously [33]. Sample DNA were amplified by PCR 
according to the manufacturer's protocol. The samples were 
genotyped using the ligation detection reaction (LDR) method 
[34] (technical support from the Biowing Applied Biotechnolo- 



gy , Shanghai, China). Analyses were repeated in 160 random 
samples (12.17%) with high DNA quality for quality control. 

Statistical Analyses 

Variations of demographic characteristics, selected variables, 
and genotypes of the OPG rs3102735 T>C, rs2073618 G>C, 
RANK rs 1805034 T>C, RANKL rs9533156 T>C and 
rs2277438 A>G variants between the cases and controls were 
evaluated using the ~/ 2 test. The associations between the five SNPs 
and risk of ESCC were assessed by computing the ORs and their 
95% CIs using logistic regression analyses for crude ORs and 
adjusted ORs when adjusting for age, sex, tobacco smoking and 
alcohol drinking status. The Hardy- Weinberg equilibrium (HWE) 
was tested by a goodness-of-fit y 2 test to compare the observed 
genotype frequencies to the expected ones among the control 
subjects. All statistical analyses were performed with SAS 9.1.3 
(SAS Institute, Cary, NC). 

Results 

Characteristics of the study population 

Characteristics of cases and controls included in the study are 
summarized in Table 1. The cases and controls are adequately 
matched on age and sex as evaluated by the y 2 tests. Notably, 
significant difference exists on both tobacco smoking and alcohol 
drinking status between the ESCC patients and the controls (p< 
0.001). 

The primary information for five genotyped SNPs was shown in 
Table 2. In all 1315 samples, the success rate of genotyping was 
95.13%, 96.35%, 96.43%, 96.43% and 96.81% for OPG 
rs3102735 C>T, OPG rs2073618 G>C, RANK rsl805034 T> 
C, RANKL rs9533156 T>C and RANKL rs2277438 A>G, 
respectively. The concordance rates of repeated analyses reached 
100%. As for the Minor Allele Frequency (MAT), there was no 
significant difference between our controls and database of 
Chinese subjects for all five SNPs. The observed genotype 
frequencies for OPG rs3102735 OT, OPG rs2073618 G>C, 
RANK rs 1805034 T>C, RANKL rs9533156 T>C and RANKL 
rs2277438 A>G polymorphisms in the controls were in HWE 
(£ = 0.191, 0.371, 0.531, 0.488 and 0.700, respectively)(Table 2). 

Associations between OPG rs31 02735 C>T, OPG 
rs2073618 G>C, RANK rs 1805034 T>C, RANKL rs9533156 
T>C, RANKL rs2277438 A>G polymorphisms and risk of 
ESCC 

Table 3 summarizes the genotype distribution of all five SNP 
polymorphisms in cases and controls. In the single locus analyses, 
the genotype frequencies of RANK rsl805034 T>C were 45.9% 
(TT), 42.9% (TC) and 1 1.2% (CC) in the case patients and 50.5% 
(TT), 41.8% (TC) and 7.7% (CC) in the control subjects. In the 
recessive model, when the RANK rs 1805034 TT/TC genotypes 
were used as the reference group, the CC homozygote genotype 
was associated with a significantly increased risk for ESCC (CC vs. 
TT/TC: adjusted OR = 1.52, 95% CI = 1.03-2.24, p = 0.036). 
When the RANK rs 1805034 TT homozygote genotype was used 
as the reference group, the TC genotype was not associated with 
the risk of ESCC (TC vs. TT: adjusted OR = 1.16, 95% CI 
= 1.03-2.24, p = 0.231), but the CC genotype was associated with 
the risk of ECSS (CC vs. TT: adjusted OR = 1.62, 95% CI 
= 1.08-2.44, p = 0.019). In the dominant model, the RANK 
rsl805034 CC variant was associated with the risk of ESCC as 
compared with the TT genotype (adjusted OR = 1.62, 95% CI 
= 1.08-2.44,^ = 0.019). 
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Table 1. Distribution of selected demographic variables and risk factors 


in ESCC cases 


and controls. 






Variable 


Cases (n 


= 629) 


Controls (n = 


686) 


ff 




n 


% 


n 


% 




Age (years) mean ± SD 


62.85 (±8 


13) 


62.58 (±7.89) 




0.541 


Age (years) 










0.155 


<63 


310 


49.28 


365 


53.21 




==63 


319 


50.72 


321 


46.79 




Sex 










0.185 


Male 


444 


70.59 


461 


67.20 




Female 


185 


29.41 


225 


32.80 




Tobacco use 










<0.001 


Never 


355 


56.44 


499 


72.74 




Ever 


274 


43.56 


187 


27.26 




Alcohol use 










<0.001 


Never 


428 


68.04 


526 


76.68 




Ever 


201 


31.96 


160 


23.32 





a Two-sided x 2 test ar, d student t test; Bold values are statistically significant (p<0.05). 
doi:1 0.1 371 /journal.pone.01 01 705.t001 



No association was detected among OPG rs3 102735 OT, 
OPG rs2073618 G>C, RANKL rs9533156 T>C, RANKL 
rs2277438 A>G polymorphisms and the risk of ECSS (Table 3). 

Stratification analyses of RANK rsl 805034 T>C genotype 
and risk of ESCC 

To evaluate the effects oiRANK rsl805034 T>C genotype on 
ESCC risk according to different age, sex, smoking and alcohol 
consumption; we performed the stratification analyses (Table 4). A 



significandy increased risk of ESCC associated with the RANK 
rsl805034 T>C polymorphism was evident among male patients 
(CCvs. TT: adjusted OR = 1.89, 95% CI = 1. 16-3.08,/) = 0.01 1) 
(TC/CC vs. TT, adjusted OR = 1.38, 95% CI =1.05-1.81, 
p = 0.022) (CC vs. TT/TC, adjusted OR = 1.68, 95% CI = 1.05- 
2.69, /? = 0.031). Likewise, in elder patients (&63 years old), 
RANK rsl 805034 T>C polymorphism was also associated with a 
significantly increased risk of ESCC (CC vs. TT, adjusted OR 
= 1.84, 95% CI = 1.02-3.31,/) = 0.041) (Table 4). 



Table 2. Primary information for OPG rs31 02735 T>C, rs2073618 G>C, RANK rsl 805034 T>C, RANKL rs9533156 T>C and 
rs2277438 A>G polymorphisms. 



Genotyped SNPs 


O/>Grs3102735 OT 


OPG rs207361 8 G>C 


RANK rsl 805034 T>C 


RANKL rs9533156 T>C 


RANKL rs2277438 A>G 


Chromosome 


8 


8 


18 


13 


13 


Gene Official Symbol 


TNFRSF1 1 B 


TNFRSF11B 


TNFRSF1 1 A 


TNFSF11 


TNFSF11 


Function 


nearGene-5 


missense 


missense 


missense 


intron region 


Chr Pos (Genome Build 36.3) 


120034251 


120033233 


58178221 


42045671 


42053168 


Regulome DB Score 3 


5 


4 


5 


5 


No Data 


TFBS b 


Y 










Splicing (ESE or ESS) 




Y 


Y 






miRNA (miRanda) 












miRNA (Sanger) 












MAF C for Chinese in database 


0.134 


0.308 


0.300 


0.439 


0.300 


MAF in our controls (n = 686) 


0.164 


0.263 


0.286 


0.464 


0.314 


p value for HWE d test 
in our controls 


0.191 


0.371 


0.531 


0.488 


0.700 


Genotyping method e 


tDR 


LDR 


tDR 


LDR 


LDR 


% Genotyping value 


95.13% 


96.35% 


96.43% 


96.43% 


96.81% 



a http://www.regulomedb.org/; 

b TFBS: Transcription Factor Binding Site (http://snpinfo.niehs.nih.gov/snpinfo/snpfunc.htm); 
C MAF: minor allele frequency, OPG rs2073618 G>C MAF is in CHB+JPT population; 
d HWE: Hardy-Weinberg equilibrium; 
e LDR: Ligation Detection Reaction. 
doi:1 0.1 371 /journal.pone.01 01 705.t002 
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Discussion 

In this hospital-based case-control study of ESCC, we investi- 
gated the association of OPG rs3102735 C>T, OPG rs2073618 
G>C, RANK rsl805034 T>C, RANKL rs9533156 T>C and 
RANKL rs2277438 A>G polymorphisms with risk of ESCC in a 
Chinese population. Our multivariable logistic analyses demon- 
strated that RANK rs 1805034 T>C genotype has an increased 
risk of ESCC. Significant association with increased risk of ESCC 
was noticed among male patients and elder patients (S63 years 
old). To our knowledge, this is the first study demonstrating a 
significant association between the RANK rsl805034 T>C 
genotype with the susceptibility of ESCC. 

OPG was initially derived from an expressed sequence tag of a 
fetal rat intestine cDNA library encoding a 40 1 -amino-acid 
polypeptide [11]. Subsequently, a physiological role of OPG in 
the maintenance of normal bone mass was underscored by several 
studies [11,35,36]. The later finding in murine myelomonocytic 
cell line 32D led to the identification of OPG binding protein or 
OPGL, which has identical sequence as RANKL and was further 
implicated with the osteoclast development [12]. Direct sequenc- 
ing of a human bone marrow-derived myeloid dendritic cell 
cDNA library identified RANK as a novel TNFR homologue 
[13]. Subsequently, RANKL was identified from murine thymoma 
cell line EL40.5 [13] as well as in T cells [37]. RANKL exists as a 
homotrimer and induces receptor clustering upon engaging 
RANK on the cell surface, consequently causes receptor 
clustering. Activation events within the cell are initiated through 
TNFR-associated factors following sufficient RANK clustering. 
Genetic variants in the OPG locus have previously been 
implicated with osteoporotic fracture [38], bone turnover [39], 
bone mineral density [40], osteonecrosis [41], diabetic neuroar- 
thropathy [42] as well as ankylosing spondylitis [43] . Alterations at 
the RANK locus and/or functionally related genes, such as 
RANKL, have also been reported to be associated with 
rheumatoid arthritis [30], aortic calcification [44], bone mineral 
density [39] and Paget's disease of bone [31]. Recently, emerging 
evidence has indicated an association between OPG/ RANK/ 
RANKL gene polymorphisms with carcinogenesis. Several studies 
demonstrated additional loci to be associated with breast cancer 
including the chromosomal region 8q24 for OPG gene [45,46]. 
SNP rs3102735 of the OPG gene has been reported to be 
associated with the susceptibility of breast cancer in Caucasian 
population [10]. Similarly, a genetic variant near the 5'-end of 
RANK (rs7226991) was associated with a breast cancer risk [47]. 
The mechanism underlying the association remains obscure so far. 
Yet, vast majority of the association on chromosome 8q24 lies at 
approximately 128 Mb and is related to various tumor entities in 
addition to breast cancer, including prostate [48] and colon cancer 
[49]. 

Among different ethnic cohorts, the frequencies of genetic 
polymorphisms vary drastically. Our study demonstrated that the 
frequency of RANK rsl805034 C was 0.286 among 686 control 
subjects in Chinese population, which is lower than that of European 
(0.438) and African American (0.478), but similar with the Japanese 
population (0.311). However, interestingly, another study reported 
the frequency of RANK rsl805034 C was 0.476 in Han population 
from North China, which differs our finding in cohort from East 
China, suggesting the ethnical impact could also be interfered with 
regional environmental factors (http://www.ncbi.nlm.nih.gov/ 
projects/SNP/snp_ref.cgi? rs= 1805034). Using Power and Sample 
Size Calculation (PS, version 3.0.43, 2009, http://biostat.mc. 
vanderbilt.edu/ twiki/bin/view/Main/PowerSampleSize), consider- 
ing RANK rs 1805034 T>C mutant alleles, the power of our analysis 
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(a = 0.05) was 0.946 in 629 ESCC cases and 686 control subjects 
with adjusted OR 1 .52. In male subgroup, the power of our analysis 
was 0.995 among 434 cases and 438 control subjects, with the 
adjusted OR 1.89. In elder cohort (—63), the power of analysis 
(a = 0.05) was 0.962 among 3 14 cases and 3 1 2 controls with adjusted 
OR 1.84. The current study has revealed the association with 
increased risk of ESCC among male patients and elder patients (&63 
years old), which was in consistent with the previous report. In a 
retrospective study involving 74,854 ESCC patients from North 
China, the prevalence among males was higher than that among 
females, similar to our findings. Moreover, this study demonstrated 
that although the prevalence significantly declined, the median age- 
of-onset of ESCC postponed [50], verifying our notion that elder 
population has higher risk. 

In conclusion, our study provides with the evidence that 
functional polymorphism of RANK rs 1805034 T>C is associated 
with the susceptibility of ESCC. We acknowledge there are several 
limitations in this study that need to be addressed. First of all, the 
study subjects were all recruited from several local medical centers 
within same area, which may not completely represent the general 
Chinese population, especially when diverse regional environmen- 
tal factors exist. Secondly, the detailed information regarding 
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